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Abstract 

Directed transport of overdamped Brownian particles driven by fractional Gaussian noises is 
investigated in asymmetrically periodic potentials. By using Langevin dynamics simulations, we 
find that rectified currents occur in the absence of any external driving forces. Unlike white 
Gaussian noises, fractional Gaussian noises can break thermodynamical equilibrium and induce 
directed transport. Remarkably, the average velocity for persistent fractional noise is opposite to 
that for anti-persistent fractional noise. The velocity increases monotonically with Hurst exponent 
for the persistent case, whereas there exists an optimal value of Hurst exponent at which the 
velocity takes its maximal value for the anti-persistent case. 
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I. INTRODUCTION 



In systems possessing spatial or dynamical symmetry breaking, Brownian motion com- 
bined with unbiased external input signals, deterministic or random alike, can assist directed 
motion of particles in nanoscale systems[l|. The subject of the fluctuation-induced transport 
was motivated by the challenge to explain unidirectional transport in biological systems {2!, 
as well as their potential technological applications ranging from classical non-equilibrium 
modelsjsl to quantum systems 4]. Ratchets have been proposed to model the unidirectional 
motion driven by zero-mean non-equilibrium fluctuations. Broadly speaking, ratchet de- 
vices fall into three categories depending on how the applied perturbation couples to the 
substrate asymmetry: rocking ratchets js], flashing ratchets jo], and correlation ratchets 3]- 
Additionally, entropic ratchets, in which Brownian particles move in a confined structure, 
instead of a potential, were also extensively studied |8|. 

Anomalous diffusion has attracted growing attention, being observed in various fields of 
physics and related sciences [9|, where by contrast with Brownian motion, long-range tem- 
poral correlations induce nonstandard dynamical behaviors. The diffusion is characterized 
through the power law form of the mean-square displacement {x'^{t)) oc t". According to 
the value of the index a, one can distinguish subdiffusion (0 < a < 1), normal diffusion 
{a = 1) and superdiffusion {a > 1). In the literature, two popular stochastic models have 
been used to account for anomalous diffusion. The first model is the continuous-time ran- 
dom walk 



14| . In this model, the subdiffusion is caused by the long waiting time between 



successive lumps and the superdiffusion is induced by the long jumps. In the minimal Levy 



ratchet 



10Nl4l|. the heavy-tailed distribution of the a-stable noise can break the thermody- 



namical equilibrium and induce directed transport. 

The second model is fractional Brownian motion (FBM) introduced by Mandelbrot and 



Van Ness 



15|. FBM has wide applications in some complex systems, such as monomer 



diffusion in a polymer chain 16|], diffusion of biopolymers in the crowded environment |l7l |. 
single file diffusion 0, and translocation of the polymer passing through a pore The 
statistical properties of FBM are characterized by the Hurst exponent < H < 1. In 
particular, its meansquared displacement satisfies (a;^(t)) oc t"^^ , thus for H < 1/2 one can 
obtain the subdiffusive dynamics, whereas for H > 1/2 the superdiffusive one 2o|. In the 



last few years, there has been a growing interest in the study of the FBM 2lN24| . However, 
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most studies of FBM focus on the free FBM and a few studies on FBM have been involved 



the potentials. Recently, Sliusarenko and coworkers 25| studied the escape from a potential 



well driven by fractional Gaussian noises and found that the escape becomes faster for de 



creasing values of Hurst exponent. Chaudhury and Cherayila [26| studied the first passage 
time distribution for barrier crossing in a double well under fractional Gaussian noises. It 
is uncertain whether fractional Gaussian noise can induce directed transport in the absence 
of any external driving forces. In order to answer this question, we studied the transport 
of overdamped Brownian particles driven by fractional Gaussian noises in asymmetrically 
periodic potentials. We focus on finding the rectified mechanism and how noise intensity 
and Hurst exponent affect the transport. 

II. MODEL AND METHODS 

In this study, we consider the directed transport of the Brownian particles driven by 
fractional Gaussian noises in the absence of whatever additional time- dependent forces. 
The overdamped dynamics can be described by the following Langevin equation in the 
dimensionless form 

dx 



V^ = -U{x) + ^VkBnH{t), (1) 

where ^nit) is the fractional Gaussian noise, D = is noise intensity, and H is the Hurst 
exponent, t] is the friction coefficient of the particle, /c^ is the Boltzmann constant, and 
T is the absolute temperature. The prime stands for differentiation over x. U{x) is an 
asymmetrically periodic potential 

U{x) = -t/o[sin(x) + ^sin(2x)], (2) 

where Uq denotes the height of the potential and A is its asymmetric parameter. 

Fractional Gaussian noise is a zero mean stationary random process with long memory 
effects 21-25|. It is closely related to the FBM process which is defined as a Gaussian 



process with an exponent < H < 1 and 



iUt)) = 0, (3) 



3 



Ai Bao-Quan, He Ya-Feng, and Zhong Wei-Rong, Phys. Rev. E 82, 061102 

(2010) 



for < s < t. In the long time limit, the autocorrelation function will decay as 

{U0)Ut))cx2Hi2H-l)t'''-\ (5) 

for < H < 1 and H ^ ^. When H = ^, fractional Gaussian noise reduces to white Gaussian 
noise. From Eq. (5), it is easy to find that the noises are positively correlated (persistent 
case) ioT ^ < H < 1, and negatively correlated (anti-persistent case) for < < |. 

Though FBM is an old topic, the consistent analytical methods are still not available. 
Here we will study the transport of the Brownian particles by using Langevin dynamics 
simulations. From Eqs. (1) and (4) one can obtain the discrete time representation of Eq. 
(1) for sufficiently small time step 6t 

x{tn+i) = x{tn) - U\x{tn))5t + VDSt^^iuin), (6) 

where n = 0A,2... and Cnin) is fractional Gaussian random number. We used the method 



described in 



25 



28 



29j for simulating fractional Gaussian random number. 
In this study, we mainly focus on the transport of the driven particles. The average 
velocity v is used to measure the transport, 

where N is the number of the realizations and to and t are the initial and the end time for 
the simulations, respectively. The asymptotic velocity V is 

V = lim v(t). (8) 



III. NUMERICAL RESULTS AND DISCUSSION 



In order to check the convergence of the algorithm, we have studied the dependence of 
average velocity on time step St, the end time t, and the number N of the realizations. From 
Fig. 1(a), (b) and (c), we can see that the algorithm is convergent and the numerical results 
do not depend on the calculation parameters {6t, t, and N) when 6t < 10^^, t > 10^, and 

> 10^. Therefore, in our simulations, the number of the realizations is more than 4 x 10'^ 
realizations, time step is chosen to be smaller than 10"'^ and t = 10^. Fig. 1(d) shows the 
estimated relative errors as a function of the noise intensity L> at = 4 x 10^, t = 10^, and 
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FIG. 1: Convergence and relative errors of the algorithm, (a) Dependence of V on time step 6t 
dX N = 10^ and t = 10^; (b)dependence of V on the end time t for simulations at St = and 
N = 10^; (c) dependence of V on the number N of the realizations at 6t = 10^^ and t = 10^; (d) 
dependence of estimated relative errors on noise intensity D. The other parameters are D = 0.5, 
Uo = 1.0, and A = 1.0. 

St = 10^'^. It is found that the relative errors are less than 0.01 even for large values of the 
noise intensity. Therefore, the parameter we used are sufficient to obtain consistent results. 

First, we study the properties of FBM for both persistent and anti-persistent cases. Figure 
2 shows the simulated sample paths for different values of H. The differences among these 
three cases are clear. For H — 0.3, the negative correlation accounts for high variability. 
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FIG. 2: (Color online) Samples of FBM for different values of the Hurst exponent H = 0.3, 0.5, 
and 0.7. 



whereas the sample is more smooth for if = 0.7 due to the positive correlation 29 1. 

The frequency spectrum of the external drive is very important to determine direction 
of motion of the Brownian particles in periodic potentials. In our minimal ratchet setup, 
fractional Gaussian noise is the only external drive, so it is necessary to analyze its frequency 
properties. In Fig. 3, we investigate the spectral density of fractional Gaussian noise for 
different values of H. We can find that the distributions of the frequency are different 
for the three cases. For white Gaussian noise {H = 0.5), the spectral density is uniform. 
However, the low frequency component is larger than the high frequency part in the spectral 
density for the persistent case {H = 0.7). For the anti-persistent case {H = 0.3), the high 
frequency component is larger in the spectral density. In order to facilitate the analysis of 
the driving mechanisms, persistent fractional Gaussian noise can be artificially divided into 
two frequency components: white Gaussian noise and low frequency ac drive. Similarly, the 
anti-persistent fractional Gaussian noise in frequency domain is equivalent to a compound 
of white Gaussian noise and high frequency drive. 

Next, we will study the directed transport mechanism for our ratchet. Usually, the 



ratchet mechanism demands three key ingredients (a) nonlinear periodic potential, 
(b)asymmetry of the potential or external driving forces, and (c)fiuctuating. Figure 4 (a) 
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FIG. 3: (Color online) The spectral density of the noises for different values of Hurst exponent 



H = 0.3, 0.5, and 0.7. 




FIG. 4: (a)Average velocity w as a function of the asymmetric parameter A of the potential for 
both persistent and anti-persistent cases at D = 0.3 and Uq = 1. (b)Asymmetrically periodic 
potential for A > and A < 0, Li is the length from the minima of the potential to the maxima 
from the left side and L2 is the length from the right side. 
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shows the average velocity as a function of the asymmetry of the potential. For the per- 
sistent case {H — 0.7), the velocity is positive for A > 0, zero at A = 0, and negative for 
A < 0. However, for anti-persistent case one can obtain the opposite velocity, negative for 
A > and positive for A < . Moreover, for both cases, there exists an optimal value of 
A at which the velocity takes its extremal value. When A — )■ 0, the system is absolutely 
symmetric and directed transport disappears. When A — )■ oo, the asymmetric potential 
described in Eq. (2) reduces to symmetric one {U{x) — — ^Asin(2x))with higher barriers, 
resulting in zero velocity. 

Now we will give the physical interpretation of the directed transport for the case of 
A = 1 (see the upper of Fig. 4 (b)). We define three time periods that are very important 
for explanation of the directed transport: driving period T, diffusion time Ti for crossing the 
steeper slope (the left side) from the minima, and diffusion time T2 for crossing the gentler 
slope(the right side). Because of Li < L2, Ti is always less than T2. Firstly, the particles 
stay in the minima of the potential (see Fig. 4 (b)) until they are catapulted out of the well 
by a large amphtude fluctuation. For the persistent case {H — 0.7), the fractional Gaussian 
noise contains more low frequency components (see the upper of the Fig. 3) and it can be 
divided into two parts: white Gaussian noise and low frequency ac drive. Due to the low 
frequency, the drive has a very long period and T ^ T2 > Ti. All particles get enough time 
to cross both sides from the minima of the potential before the drive reverses its direction. 
However, the left side is steeper than the right one, more particles climb the barrier from 
the right side, so the average velocity is positive. For the anti-persistent case{H — 0.3), 
the high frequency components dominate over the low frequency ones(see the bottom of the 
Fig. 3). In this case, the drive has a short period and Ti < T < T2 (or T < Ti < T2). In 
a short driving period, the particles have sufficient time to diffuse across the steeper side of 
the well, resulting in negative average velocity. It should be pointed out that the average 
velocity will tend to zero for the case of T -C Ti <T2 (very small values of H) which is also 
shown in Fig. 6. 

The noise intensity dependence of the average velocity is shown in Fig. 5 for different 
values of Hurst exponent. The curve is observed to be bell shaped. When D — > 0, the 

particles cannot pass across the barrier and there is no directed current. When Z) — )■ 00 so 
that the noise is very large, the effect of the potential disappears and the average velocity 
tends to zero, also. Therefore, one can see that the curves demonstrate nonmonotonic 
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FIG. 5: Average velocity as a function of noise intensity D for different values of Hurst exponent 
H = 0.2, 0.3, 0.6, 0.7, and 0.8 at A = 1.0 and Uq = 1. 



behavior. 

Figure 6 displays the Hurst exponent dependence of the average velocity at A = 1.0. It 
is found that the average velocity is positive ioi H > |, zero at H = |, and negative for 
H < ^. For the persistent case {H > |), the average velocity increases monotonically with 
the Hurst exponent. However, for the anti-persistent case {H < |), there exists a value of 
Hurst exponent at which the average velocity takes its extremal value. When if — )■ 0, the 
noise term in Eq. (1) will disappear, the system is deterministic and the directed transport 
also disappears. When — > |, the fractional Gaussian noise reduces to the white Gaussian 
noise, the system undergoes thermal equilibrium and the average velocity tends to zero. 



IV. CONCLUDING REMARKS 



In this paper, we studied the directed transport of overdamped Brownian particles driven 
by fractional Gaussian noises. From numerical simulations, we can find that fractional Gaus- 
sian noise can break the detailed balance and induce directed transport. Similar to the 
classic ratchets ij , there exists a value of noise intensity at which the average velocity takes 



its extremal value. The average velocity as a function of the asymmetry of the potential is 
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FIG. 6: Average velocity V as a function of the Hurst exponent H at A = 1.0, D = 0.5, and 
Uo = l. 



monotonic. From the numerical analysis of the spectral density of fractional Gaussian noises, 
it is found that the low frequency component in spectral density is larger than the high fre- 
quency component for the persistent case {H > |), whereas the high frequency component is 
dominated for the anti-persistent case {H < i). Due to the difference of the spectral density 
between the persistent and anti-persistent cases, the average velocity has the opposite sign 
for the two cases. Remarkably, the average velocity increases monotonically with the Hurst 
exponent for the persistent case. However, for anti-persistent case, there exists an optimal 
value of the Hurst exponent at which the velocity takes its extremal value. 

Directed transport in static ratchet potentials can also be induced by the other types of 
noise, such as a-stable noise (Levy ratchet) ll|- IJ], white shot noise (shot-noise ratchet) 30| . 



In the Levy ratchet 



ll|-|l4|. the thermo- 



and two correlated noises (correlated ratchet) 
dynamical equilibrium is broken by the the heavy-tailed distribution of the a-stable noise. 
For shot-noise ratchet 30[, the temporal asymmetry of white shot noise can induce an ef- 
fective, inhomogeneous diffusion, so the net current occurs. In the correlated ratchet [t], 
fluctuation-induced transport is driven by both additive Gaussian white noise and additive 
colored noise. The additive colored noise can be treated as the multiplicative noise by in- 
troducing a new auxiliary variable, therefore an effective, inhomogeneous diffusion appears. 
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However, in our fractional Gaussian noise-induced ratchet, the directed transport is induced 
by the asymmetry of noise spectral density When if > |, the fractional Gaussian noise con- 
tains more low frequency components, whereas the high frequency component is dominated 
for < |. 
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